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Masses for {Qq}{Q'^ molecular states are systematically studied in QCD sum rules. The 
interpolating currents representing the related molecular states are proposed. Technically, contribu- 
tions of the operators up to dimension six are included in operator product expansion (OPE). Mass 
spectra for molecular states with {Qq}{Q^ ^q} configurations are obtained. 
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I. INTRODUCTION 

The field of heavy hadron spectroscopy has attracted much attention nowadays. Experimentally, plentiful 
hadronic resonances have been observed, such as the so-called X, Y, and Z states 

(for experimental reviews, e.g., see ^ [IB])- In theory, different from conventional charmonium states, 
some of these resonances may not reconcile with the quark model picture, hence it is not easy to find 
appropriate positions for them in mesonic spectroscopy. Many authors [ill . [Tl . Hi, Q E [li tend to 
interpret some of the hadrons as possible mesonic molecule candidates for their masses are very close to 
the meson-meson thresholds. For example, it has been proposed that X{3872) could be a D*D molecular 
state [ijj, Z+(4430) be a Li*i)imolecular state [li], ^(3930) be a D*D* molecular state [lllli, ^(4140) 
be a D*D* molecular state 15|, |l6| etc.. As a matter of fact, the charmed molecular states were put forward 
long ago in [l^ , and it has also predicted that the molecular states involving hidden cc pair do exist and 
have a rich spectroscopy in 18|. What is also very important, the existence of molecular states is not 
excluded by QCD itself. Indubitably, investigations of molecular states could deepen one's understanding 
of quark-gluon interaction. Thus, the quantitative descriptions of molecular states' properties like masses 
are quite needed for well comprehending their structures. In a word, it is interesting and significative to 
study mass spectra for the molecular states. 

Motivated by the above reasons, we devote to obtain the spectra for molecular states with {Qq}{Q^ 'g} 
configurations in this work. However, it is a great challenge to extract information on the spectrum from 
first principles. While QCD has been widely accepted as the correct theory describing strong interaction, it 
is still far from clear how to gain hadron masses from the simple Lagrangian. That's because QCD is highly 
noperturbative in the low energy region where futile to attempt perturbative calculations. Therefore, one 
has to deal with a genuinely strong field in nonperturbative methods. Before thoroughly grasping the 
essence of QCD confinement, for the moment, one could make use of QCD sum rule (19l] (for reviews see 
""3, 21, 22, 3] and references therein), which is a comprehensive and reliable working tool for evaluating 



the nonperturbative effects. The basic idea of this approach is bridging the gap between the perturbative 
and nonperturbative sectors by employing the language of dispersion relations. Already, there have been 



several works testing some molecular states from QCD sum rules up to now [2J, |25|, l26| . At present, we'd 
like to gain the spectra for \} molecular states with QCD sum rules, which could serve as a 

extension of our previous work on {Qs}{Q^ ^s} states ^22\. 

The paper's framework is as follows. In Sec. HH QCD sum rules for the molecular states are introduced, 
and both the phenomenological representation and QCD side are derived, followed by the numerical analysis 
to extract the hadronic masses and a brief summary in Sec. IIIII 



II. MOLECULAR STATE QCD SUM RULES 

The QCD sum rule attempts to link the hadron phenomenology with the interactions of quarks and 
gluons, the elementary point of which is the choice of interpolating current. Following the standard scheme 
[lC|, the Qq mesons with = 0^, 1^, 0+, and 1+ are named D, D* , Dq, and Di for charmed mesons, 
with B, B* , Bq, and Bi for bottom mesons, respectively. Here the configurations for these mesons are 
represented as (Qq), (Qq)*, (Qq)n, and {Qq)i. In full theory, the interpolating currents for these mesons 
can be found in Refs. [27|, |28|. Consequently, the interpolating currents for the related molecular states 
are constructed in following forms, with 

HQq)iQ'q) = {^ailbQa){Q' bilbqh), 

HQq)'{Q'q)' = iqal^,Qa){Q'b7''qb), 
MQq)UQ'q)o = i^aQa){Q'b<lb), 
MQq)l{Q'q)l = iQalt,l5Qa)iQ'bl^l5qb), 
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for one type of hadrons, and 



i(Q9)(Q'9)S = {<lah5Qa){Q'bqb), 



■^'W3)iW'«)5 = (9a7^75Qa)(Q'bgb), 



for another type, where q indicates the light quark u or d, Q and Q' denote heavy quarks {Q = Q' or 
Q 7^ Q')^ with a and b are color indices. 

For the former case, the starting point is the two-point correlator 

U{q^) = i J d^xe"i-^0\T\j{x)j+m\0). (1) 

The correlator can be phenomenologically expressed as a dispersion integral over a physical spectral func- 
tion 



Mlj-q2 TT J,^ s - q 



n(g^) = " + - ds + subtractions, (2) 



where Mh denotes the mass of the hadronic resonance, and A^f gives the coupling of the current to the 
hadron {Q\j\H) = Xh- In the OPE side, the correlator can be written in terms of a dispersion relation as 

n(g2) = / ds ^ {mQ = mQ, or mQ^niQ,), (3) 

where the spectral density is given by the imaginary part of the correlator 

TT 

After equating the two sides, assuming quark-hadron duality, and making a Borel transform, the sum rule 
can be written as 

y2^^-Ml/M^ = dsp'^^^{s)e-^'^\ (5) 

where indicates Borel parameter. To eliminate the hadronic coupling constant Xh, one reckons the 
ratio of derivative of the sum rule and itself, and then yields 



For the latter case, one starts from the two-point correlator 

W^{q^) = I I d\e''i-^0\T[3^{x)f+{0m. (7) 



Lorentz covariance implies that the two-point correlation function can be generally parameterized as 



n^-(g2) = _ g>^-)Ii^^){q^) + ^nW(g2). (8) 
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The part of the correlator proportional to g^i, will be chosen to extract the mass sum rule here. In 
phenomenology, n(^^(g^) can be expressed as a dispersion integral over a physical spectral function 



n^^Uq^) = \ ' „ + - / ds 5— + subtractions, (9) 

-q^ TT J,^ s- q^ 



In the OPE side, H^^^q"^) can be written in terms of a dispersion relation as 

^OPE^ 
ds- 

l{mQ+mQ, 

where the spectral density is given by 

pOPE(^) ^ (11) 



Tl^^\q^)= dsf^ Y {niQ^mQ, oimQ^mq,), (10) 



TT 



After equating the two sides, assuming quark-hadron duality, and making a Borel transform, the sum rule 
can be written as 

^yi^i)^2^-Ml/M- = r dsp^'^^{s)e-'"'\ (12) 

J {niQ+mQ,)'^ 

To eliminate the hadronic coupling constant \^^\ one reckons the ratio of derivative of the sum rule and 
itself, and then yields 

Ml ^ r dsp^^^se-^'^'" I r dspOP\-^/M\ (13) 

Calculating the OPE side, one works at leadirig order in and considers condensates up to dimen- 
sion six with the similar techniques in Refs. 2^, 3^. To keep the heavy-quark mass finite, one uses the 
momentum-space expression for the heavy-quark propagator. One calculates the light-quark part of the 
correlation function in the coordinate space, which is then Fourier-transformed to the momentum space 
in D dimension. The resulting light-quark part is combined with the heavy-quark part before it is di- 
mensionally regularized at Z? = 4. For the heavy-quark propagator with two and three gluons attached, 
the momentum-space expressions given in Ref. 27| are used. After some lengthy OPE calculations, the 
concrete spectral densities can be acquired, which are collected in the Appendix. 



III. NUMERICAL ANALYSIS 

In this section, the sum rules (El) and (|13p will be numerically analyzed. The input values are taken as 
rric = 1.23 GeV, mt = 4.20 GeV with (qq) = -(0.23)3 GeV^ {gqa ■ Gq) = ml {qq), ml = 0.8 GeV^ 
(g2G2) ^ 0.88 GeV, and {g^G^) = 0.045 GeV^ 00,0. Complying with the standard procedure of 



sum rule analysis, the threshold sq and Borel parameter are varied to find the optimal stability window, 
in which the perturbative contribution should be larger than the condensate contributions while the pole 
contribution larger than continuum contribution. Thus, the regions of thresholds are taken as those values 
presented in the related figure captions, with AP = 3.5 - 4.5 GeV^ for DD, D*D, D*D*, D^Dl, DiD^, 
DD^, DiD, D*Dl, and D*Di, = 7.5 - 9.0 GeV^ for BD, B*D, B*D*, B^D^, BiD^, BiDi, 
BD*, BiD, B*D*, B*Di, D*B, DiB*, DB*, DiB, D*B*, and D*Bi, and = 9.5 - 11.0 GeV^ for 
BB, B*B, B*B*, B^B^, BiB^, BiBi, BB^, BiB, B*B^, and B*Bi, respectively. The corresponding 
Borel curves are exhibited in Figs. 1-18, and the numerical results are listed in Tables I-II. It is worth 
noting that the numerical errors reflect the uncertainty due to sum rule windows (variation of the threshold 
So and Borel parameter M^) only; the uncertainty resulted from the variation of the quark masses and 
QCD parameters is not included. After a comparison, one can find that the numerical results for D*D, 
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TABLE I: The mass spectra of molecular states with same heavy quarks. 



Hadron 


configuration 


mass (GeV) 


Hadron 


configuration 


mass (GeV) 


DD 


{cq){cq) 


3.76±0.f0 


BB 


{bq){bq) 


f0.58±0.f0 


D'D 


{cqYicq) 


3.88±0.f0 


B'B 


{bqy(bq) 


f0.62±0.f0 


0*0* 


{cq)*{cq)* 


3.9f ± O.ff 


B*B'' 


(6g)*(feg)* 


f0.67±0.f0 


Do Do 


(cq)o(cq)o 


4.56±0.ff 




{bq)o(bq)o 


ff.28±0.08 


DiDo 


(cq)i(cq)n 

\ ^ j ^ \ ^ ) u 


4.62 ± O.ff 


BiBq 


(bq)i(bq)n 


ff.32±0.09 


DiDi 


(cg)i(cg)i 


4.66±0.f3 


BiBi 


{bq)i(bq)i 


ff.33±0.f2 


DDI, 


(cg)(cg)S 


4.2f ±0.07 


bb:, 


{bq)(bq)o 


ff.03±0.09 


DiD 


(cg)i(cg) 


4.34 ± 0.07 


BiB 


ibqUh) 


ff.04±0.09 




(eg)* (eg),* 


4.26 ±0.07 


B*B*o 


{bqy{bq)o 


ff.02±0.09 


D*Di 


(cg)*(cg)i 


4.44 ± 0.09 


B*Bi 


{bqnbq)i 


ff.03±0.09 



TABLE ff: The mass spectra of molecular states with differently heavy quarks. 



Hadron 


configuration 


mass (GeV) 


Hadron 


configuration 


mass (GeV) 


BD 


{bq){cq) 


7.f2±0.09 


B*D*o 


(6g)*(cg)(* 


7.67 ±0.06 


B*D 


{bq)*{cq) 


7.28 ±0.09 


B*Di 


(6g)*(cg)i 


7.74 ±0.07 


B'D* 


{bq)*{cqy 


7.29±0.f0 


D*B 


(cg)*(6g) 


7.2f ±0.09 


B*oD'o 


(6g)o(cg)5 


8.04 ±0.08 


DiB:, 


(eg)i(&g)S 


8.04±0.f0 


BiA) 


(6g)i(cg)5 


8.06±0.f3 


DBl, 


{cq)(bq)o 


7.70 ± 0.06 


BiDi 


(6g)i(cg)i 


8.07 ± O.ff 


DiB 


(cg)i(6g) 


7.74 ±0.07 


BD'o 


{bq){cq)o 


7.68 ± 0.06 


D*B*o 


(cg)*(6g)(* 


7.76 ± 0.06 


BiD 


(6g)i(cg) 


7.77 ±0.06 


D'Bi 


(cg)*(6g)i 


7.76 ± 0.07 



D*Di, and D*D* are in good agreement with the experimental data for X(3872), Z+(4430), and y(3930), 
respectively, which could support the molecular interpretations for these hadrons. 




FIG. f : The dependence on for the masses of DD and BB from sum rule ([6|. The continuum thresholds are 
taken as — 4.f ~ 4.3 GeV and = fO.9 ~ f f .f GeV, respectively. 



6 




FIG. 2: The dependence on for the masses of D* D and B*B from sum rule (|13p. The continuum thresholds 
are taken as ^/so = 4.3 ~ 4.5 GeV and ^/so = 11.0 ~ 11.2 GeV, respectively. 




FIG. 3: The dependence on for the masses of D* D* and B* B* from sum rule ([6]). The continuum thresholds 
are taken as ^/sq = 4.3 ~ 4.5 GeV and ^/so = 11.0 ~ 11.2 GeV, respectively. 




FIG. 4: The dependence on for the masses of DiDq and BiBq from sum rule (|13p . The continuum thresholds 
are taken as = 5.0 ~ 5.2 GeV and y/so = 11.6 ~ 11.8 GeV, respectively. 
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FIG. 5: The dependence on for the masses of DqDq and BqBq from sum rule ([6]). The continuum thresholds 
are taken as ^/so = 4.9 ~ 5.1 GeV and ^/so = 11.5 ~ 11.7 GeV, respectively. 
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-5.4GeV 
-5.5GeV 
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FIG. 6: The dependence on for the masses of Di-Di and BiB\ from sum rule ([6]). The continuum thresholds 
are taken as ^/sq = 5.3 ~ 5.5 GeV and = 12.0 ~ 12.2 GeV, respectively. 
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FIG. 7: The dependence on for the masses of DD^ and BBq from sum rule (O. The continuum thresholds 
are taken as ^/sq = 4.5 ~ 4.7 GeV and = 11.4 ~ 11.6 GeV, respectively. 
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FIG. 8: The dependence on for the masses of DiD and BiB from sum rule (|13p . The continuum thresholds 
are taken as ^/so = 4.7 ~ 4.9 GeV and ^/so = 11.4 ~ 11.6 GeV, respectively. 
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FIG. 9: The dependence on for the masses of D*Dq and B* Bq from sum rule H13|l . The continuum thresholds 
are taken as = 4.6 ~ 4.8 GeV and ^/so = 11.4 ~ 11.6 GeV, respectively. 
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FIG. 10: The dependence on for the masses of D* D\ and B*Bi from sum rule ([6]). The continuum thresholds 
are taken as = 4.9 ~ 5.1 GeV and y/so = 11.4 ~ 11.6 GeV, respectively. 
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FIG. 11: The dependence on for the masses of BD and B*D* from sum rule (|6]). The continuum thresholds 
are taken as = 7.5 ~ 7.7 GeV and ^/so = 7.7 ~ 7.9 GeV, respectively. 
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FIG. 12: The dependence on for the masses of BqDq and B\D\ from sum rule (|6]). The continuum thresholds 
are taken as ^/sq = 8.3 ~ 8.5 GeV and = 8.6 ~ 8.8 GeV, respectively. 
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FIG. 13: The dependence on for the masses of D* B and B*D from sum rule (|13p . The continuum thresholds 
are taken as ^/so = 7.6 ~ 7.8 GeV and = 7.7 ~ 7.9 GeV, respectively. 
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FIG. 14: The dependence on for the masses of D\Bq and B\Dq from sum rule (|13p . The continuum thresholds 
are taken as = 8.4 ~ 8.6 GeV and y/ao = 8.4 ~ 8.6 GeV, respectively. 
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FIG. 15: The dependence on for the masses of D*Bq and B*Dq from sum rule (I13|) . The continuum thresholds 
are taken as = 8.1 ~ 8.3 GeV and = 8.0 ~ 8.2 GeV, respectively. 



8.1 GeV 
-8.2GeV 
-8.3GeV 



7.5 8 M=(GeV')8.5 



8.1 GeV 
-8.2GeV 
-8.3GeV 



7.5 8 M=(GeV") 8.5 



FIG. 16: The dependence on for the masses of D* Bi and B* D\ from sum rule ([6]). The continuum thresholds 
are taken as ^/so = 8.1 ~ 8.3 GeV and = 8.1 ~ 8.3 GeV, respectively. 
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FIG. 17: The dependence on for the masses of DBq and BDq from sum rule (|6]). The continuum thresholds 
are taken as = 8.0 ~ 8.2 GeV and y/ao = 8.0 ~ 8.2 GeV, respectively. 
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FIG. 18: The dependence on for the masses of D\B and BiD from sum rule (|13p . The continuum thresholds 
are taken as ^/sq = 8.1 ~ 8.3 GeV and y/so = 8.1 ~ 8.3 GeV, respectively. 

In summary, the QCD smn rules have been employed to calculate the masses for {Qq\{Q^ molecular 
states, including the contributions of the operators up to dimension six in OPE. We have attained mass 
spectra for {Qq}{Q'^ \} molecular states in the end. In molecular pictures, the numerical results for the 
masses of X(3872), Z+(4430), and F(3930) agree well with their corresponding experimental values, which 
can support that X(3872) could be a D*D molecular state, Z+(4430) be a D*Di molecular state, and 
y(3930) be a D*D* one. On all accounts, all the numerical results are expecting further experimental 
identification and it is looking forward to more experimental evidence on molecular states. 



APPENDIX 



The spectral densities are distinguished for two kinds of doubly heavy molecular states, namely, con- 
taining the same or differently heavy quarks. It is defined that r{mQ^,mQ^) — ciitiq_^ + Pttiq^ — a(3s. 
Concretely, r{mQ, mq) = aroq + Pttiq — a^s and r{mQ, mqi) — arriQ + Pttiq, — a/Ss. First, with 



pOPE(^) ^ ^pert(^) _^^(,,)(^) _^^(<,,..G,)(^) _^.^(9=G=)(^) ^^te^G^)(^)^ 



(aa)^ s ^im) r-^-^ da f^-°' dp . ,2 

P^'^^'^Hs) = £"^^^^J^~\p{l-a-P)r{mQ,mQ), 

Otrnin *' Pmin 

for (Qg)(Q(z), 

pOPE(^^ ^ _{pPert(g^ ^^(,-,)(^) ^^{9g>^(^) ^^(gg..G9>(^) +^(g^G=>(g) + ^(/G^)(g^ 

CXrnin *' Pmin '"^ 

tirriin Prmn 

P'^'^^i^) = ^ dP{l-a-(3){l+a + (3)[r{mQ,mQ) + 2ml(i], 
for (Qg)*(Q(z), 
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X /•"— da [^~" dp . ,2 



d(3{l-a- P)[r{mQ,mQ) + 2m^/3], 
for (Qg) 



P'"''''\^) = / ^ / d0{l-a-(3MmQ,mQ)+2ml/3], 

V Oirn.i-n. IJm.in. 



2l27r6 

for 



Dert/ N 3 r'-'^'^ da f^-" dp , ^, 



m( ^ da df3 ^ 



for (Qg)i(09)i, 

pOPE(^) ^ _{pPert(^) +^(,-,>(^) +^(,,>^(^) +p(s«-a.G,>(^) +^(/G^>(^) +^(s3G3>(^)^^ 



,„-„^ , , 3(gg) r-^-'' da f^"^ d(3 ^, , .2 



p(9'G^){s) = -^J^m^Q T""' ^ C "d/?(l-a-/3)(l + a + /3)r(mQ,mQ), 
for {Qq)i{Qq)l, 

P^^'^^'H^) = £'"°^5/^''"d/3(l-a-/3)r(mg,mQ), 
/^^'''''^(«) = P"'5 r\(i{l-a-P)[r{mQ,mQ)+2ml(3\, 

for {Qq){Qqro, 

pOPE(s) = _{^pert(^) ^^(,,)(^) ^^(g5..G,>(^) ^^(g^G=>(^) ^^(g3G3)(^)^^ 

p(.^G=>(,) = _.^!)^2^£"°^ "d/3(l-a-/3)(l + a + /J)r(mQ,mQ), 



for {Qq)i{Qq), 

pOPE(^) ^ _|^pert(^) +p(,,>(^) +^(«-«>^(^) +p(.,..G,>(^) +^<,=G^>(^) +p(.3G3> 
/'P'^'(^) = / 5 / i(l-«-/3)(l + " + /3MmQ,mc3)^ 

P^b''G-){s) = -^i^ r "(i/3(l-a-/3)(l + a + /3)[r(mQ,mQ) + 2m^/3], 

for {Qq)*{Qq)o, and 
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for {Qq)*{Qq)i. The integration limits are given by amin = (1 - yi^-4rn^/s)/2, Umax = (1 + 

^1- 4m^/s)/2, and = am^/isa -rn^). 
Second, with 



Dert. N 3 f"'^"" da f^'" dp ,^ 



277r4 



V OCrnAn. IJm.in. ' 



^^^^^H^) = day^^^^d/^(l-a-/3)[(^ + ^MmQ,mQ0+2(^ + -|-)], 

for (Qg)(Q'g), 



pOPE(^) ^ _{^pert(^) +^(,,)(^) +^(,,)^^^^ +^(,,-..G,>(^) +^(/G^>(^) +^(,3G3>(^)^ 



P<m)'(s) = -^|2-mQmQ/^(s-m^ + m^,)2 -4m^,s/s, 



= ^^^1^ jy^'daxi-'^ j"~"dpr{mQ,mQ,) + {^ + ^)[aml + {l-a)ml,-a{l-a)s]}, 
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p(3'G^){s) = --^^ P^'da r "d/3(l-a-/3)(l + a + /3)(!^ + !j)r(mQ,mQO, 

P^^'^n^) = da d0{l-a-P){l + a + P)[{- + -)r{mQ,mQ,) + 2i^ + ^)], 

for (QqriQ'q), 

pOPE(s) = pPert(^^ ^^(,-,>(^) ^^(«-«)^(^^ ^^<g,..G«)(^) ^^(S=G=>(^) ^^(S^G^)(^)^ 



3(g(z) 



267r4 



/ 3^3\ (q^G^) f""""" r^~" 1 1 rn%,p m%a 

P^'^'^i') = ^ da d0{l-a-m^+gs)r{mQ,mQ,)+2{^ + ^)], 

for {QqnQ'q)*, 



Dert/ N 3 r"'"^ da f^'" d(3 



(Xmin Pmin ^ ^ 



" OCmin 

for (Qg)5(Q'g)S, 

pOPE(^^ ^ pPCrt(^) ^^(«,)(^^ ^^{,-g>^(^) ^^to..G,>(^^ ^^(g^G=>(g) ^^(S^G^>(^^^ 



/„2/^2\ foc^a.^c pl — a ^2 ^2 



for (Qg)i(Q'«)i, 



OPE 



(S) = -{pPert(,) +^(,-,)(^) +^(,-,>^(,) +^(fl</<x Gg>(^) ^^(g^G=)^^^ ^^(S^G3)(^)|^ 
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212^6/ '^«/ -a -/?)(!+ a + /3)(^ + ^MmQ,mQ0, 



/ da d/3(l-a-/3)(l + a + /3)[(_+ )r(mQ,mQ0+2(- 

v OCrn.in. IJm.in. ' 



tor {QtMQ'lK, 

I- \2 

*' Otrnin *' Pmin ^ 



^^'^''H^) = da dP{l-a-m^+g-,)r{mQ,mQ>)+2{^ + ^)], 

for (Qg)(Q'g)S, 

pOPE(^) ^ _{ppert(^) +^(,-,)(^) +^(,,>=(^) +^(.,-..G,>(^) +^(,^G^>(^) +^(^^ 
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P^'^'Hs) = / da dP{l-a~P)il+a + P){^ + ^)r{mQ,mQ,), 

pio^G^)^s) = r^'^daf "d/3(l-a-/3)(l+a + /3)[(^ + ^MmQ,mQ0+2(^ + ^)], 

for {Qq)i{Q'q), 

pOPE(,) = _{pPert(^) +^(,-«>(^) +^(,-«>=(^) +^(s,-<..G,>(^) +^(,^G^>(^) +^(,3G3>(^^ 



2%2 



'^^'''''Hs) = r^^daf "rf/3(l-a-/3)(l + a + /3)(^ + ^MmQ,mQO, 

P^'^'Ks) = / da dP{l-a~P){l+a + P)[{^ + -^)r{mQ,mQ,) + 2{^ + ^)], 

for (Qg)*(Q'g)o, and 

pO^^is) = pPert(,) +^(?.>(5) +^(95-G,>(^) ^^(g^G=>(^) ^^(/G3>(^^^ 



22 



P^'''^^) = rfay^_d/3(l-a-/3)[(-3+^MmQ,mQ,) + 2(^ + -|-)], 

for (Qq)* {Q'q)i- The integration limits are given by amin = ~ "^q + ™q' ~ 

^ (s - TOq + m|,)2 - 4m|,s]/(2s), = [s - m| + + y^(s - +TOq/)^ - 4m^,s]/(2s), and 
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